IGC-11/1-2 



Loop Quantum Gravity and Cosmology: 
A dynamical introduction 

Martin BojowalcE] 

Institute for Gravitation and the Cosmos, 
The Pennsylvania State University, 
104 Davey Lab, University Park, PA 16802, USA 

Abstract 

Loop quantum gravity and cosmology are reviewed with an emphasis on evalu- 
ating the dynamics, rather than constructing it. The three crucial parts of such an 
analysis are (i) deriving effective equations, (ii) controlling the theory's microscopic 
degrees of freedom that lead to its spatial discreteness and refinement, and (hi) en- 
suring consistency and anomaly- freedom. All three issues are crucial for making the 
theory testable by conceptual and observational means, and they remain challenging. 
Throughout this review, the Hamiltonian nature of the theory will play a large role 
for properties of space-time structure within the framework discussed. 



u It would be permissible to look upon the Hamiltonian form as the fundamental 
one, and there would then be no fundamental four- dimensional symmetry in the 
theory. One would have a Hamiltonian built up from four weekly [sic] vanish- 
ing functions, given by [the Hamiltonian and diffeomorphism constraints] . The 
usual requirement of four- dimensional symmetry in physical laws would then 
get replaced by the requirement that the functions have weakly vanishing P.b. 's, 
so that they can be provided with arbitrary coefficients in the equations of mo- 
tion, corresponding to an arbitrary motion of the surface on which the state is 
defined." 
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1 Introduction 



In its different incarnations, quantum gravity must face a diverse set of fascinating problems 
and difficulties, a set of issues best seen as both challenges and opportunities. One of 
the main problems in canonical approaches, for instance, is the issue of anomalies in the 
gauge algebra underlying space-time covariance. Classically, the gauge generators, given 
by constraints, have weakly vanishing Poisson brackets with one another: they vanish 
when the constraints are satisfied. After quantization, the same behavior must be realized 
for commutators of quantum constraints (or for Poisson brackets of effective constraints), 
or else the theory becomes inconsistent due to gauge anomalies. If and how canonical 
quantum gravity can be obtained in an anomaly-free way is an important question, not yet 
convincingly addressed in full generality. Posing one of the main obstacles to a complete 
formulation of quantum gravity, this issue is hindering progress toward a detailed evaluation 
of quantum gravitational dynamics. A reliable phenomeno logical analysis must, after all, 
start with a consistent set of sufficiently general dynamical equations. 

But the strong and tough requirement of anomaly-freedom is also an opportunity, for it 
allows an analysis of quantum space-time and the changes in its structure possibly implied 
by quantum gravity. Addressing the anomaly problem is, moreover, crucial for an under- 
standing of the dynamics of quantum gravity, both in the sense of constructing consistent 
dynamical equations at the quantum level and in the sense of evaluating equations and 
their solutions to bring out physical effects. 

Although the anomaly problem has not been addressed in full generality, sev eral model 
system s have by now been analyzed in loop quantum cosmology, as reviewed by iBojowald 
(|2008cJ), with this question in mind. Loop quantum cosmology is a rather wide area within 
loop quantum gravity, analyzing several classe s of mod e l syst e ms and perturbations around 



them. Lo op quantum grav ity, detailed by iRovellil ( 120041 ). lAshtekar and Lewandowski 
(120041 ) and iThiemannl (120071 ). is a canonical quantization of general relativity based on 
holonomies (the eponymous loops) as elementary variables. The use of holonomies al- 
lows a background-independent formulation free of auxiliary metrics, and it implies several 
specific properties of the resulting dynamics. 

In all the systems used in loop quantum cosmology, quantization techniques close 
to those of a general loop quantization are used; they can thus be seen as capturing 
at least some of the crucial properties of full loop quantum gravity. To different de- 
grees, most of these models m ake additional use of symmetry reduction as introduced by 



Bojowald and Kastrupl (120001 ). simplifying much of the quantum geometry and thereby 



providing rather direct access to the much less understood quantum dynamics. Thanks to 
these steps, impli cations of the quan t um dy nam ics of loop quantum gravity, as generally 
defined based on iRovelli and Smolinl ( 119941 ) and IThiemannl ( 1l998af ) , have been evaluated 
quite explicitly for the first time. 

In general terms as well as for particular questions arising in loop quantum cosmology 
and loop quantum gravity, three key issues regarding quantum space-time arise, not unlike 
what one would expect for any fundamental gauge theory with microscopic degrees of 
freedom: 
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Effective dynamics: In quantum gravity, geometry is described unsharply by whole 
states with all their fluctuations and correlations, in addition to the expectation 
values for an average geometry Equations of motion for expectation values receive 
quantum corrections in their effective dynamics, as it may describe a quantum ge- 
ometry. Along with this consequence of quantizing gravity come not only new math- 
ematical space-time structures but also a vast enlargement of the number of degrees 
of freedom by quantum variables. 

The strongest control of such a high-dimensional dynamical quantum system is usu- 
ally obtained for dynamical coherent states, defined as states saturating uncertainty 
relations at all times. If such states exist, they provide insights into the minimal 
deviations from classical behavior expected for a quantum system. The form and be- 
havior of dynamical coherent states in loop quantum cosmology can be highlighted 
in several models, bringing out the role of space-time fluctuations and correlations 
as degrees of freedom beyond the classical ones. Exact dynamical coherent states 
exist only in special models and for specific initial values. Nevertheless, they allow 
interesting views on the generic quantum dynamics as it arises in quantum gravity. 

Before all corrections are derived for a large class of models, a clear analysis provided 
by dynamical coherent states, when they exist, unambiguously shows the first devia- 
tions from classical behavior. More generally, when exact dynamical coherent states 
do not exist, additional quantum corrections will result. They can often be computed 
perturbatively, analogously to loop corrections in interacting quantum field theories. 

Discrete dynamics: In addition to those generic effects due to non-classical state pa- 
rameters, underlying space-time structures typically change even for the expectation 
values of a quantum-gravity state. Most importantly, discrete geometry, at least 
in purely spatial terms whic h is b y now well understood in loop quantum gravity 



in purely spatial terms wnic n is b y now well understood m loop quantum gravity 
following Rovelli and Smolin ( 1995 ); Ashtekar and Lewandowski (Il997l . 1998 ). shows 



several detailed properties of importance for the dynamics and thus for space-time 
geometry. 

A spatial slice in space-time is equipped with a discrete quantum geometry, roughly 
seen as making space built from atomic patches of certain discrete sizes. One of the 
main problems of quantum dynamics is to show how these spatial atoms along slices 
fit together to form a quantum space-time — or, more dynamically, how the spatial 
atoms change, merge, subdivide and interact as time is let loose. For an expanding 
universe, one would expect the discrete spatial structure to be refined as the volume 
increases; otherwise discrete sizes would be enlarged by huge factors, especially during 
inflation, making them macroscopic. The full dynamics of loop quantum gravity 
has indeed provided several hints that the number of discrete building blocks must 
change from slice to slice, once the dynamics is implemented consistently. This lattice 
refinement can be modelled even in the simplest, most highly s ymmetric situations 



of loop quantum cosmology, laid out by iBojowaldl (120061 . l2008al ) . And it has shown 



several specific implications by which its precise form can already be constrained. 
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Consistent dynamics: Dynamics unfolds in time, but time is relative. Making sure 
that descriptions using different notions of time, corresponding to measurements by 
different observers, can agree about their physical insights requires the consistency 
conditions of general covariance. Since general covariance in gravity is implemented 
by gauge transformations, any quantization or even just a modification of the theory 
must, for consistency, respect this principle and be anomaly free. 

While the previous two points manifest themselves already in homogeneous models, 
where they can most easily be studied, the anomaly problem arises only in inhomo- 
geneous situations. Within homogeneity, all but one of the gauge transformations 
underlying covariance are fixed. Anomalies can only arise if there are at least two 
independent gauge transformations; after quantization, they would be anomalous if 
their composition is no longer a gauge transformation. One could, of course, make 
it a gauge transformation by definition, by declaring the whole group generated by 
the independent gauge transformations as the gauge group. But if this group is too 
large, it would identify variables which are to be considered physically distinct, re- 
moving observables and degrees of freedom and in many cases leaving no non-trivial 
solutions. Changing the gauge transformations of a classical theory by quantum ef- 
fects requires much care; only so-called consistent deformations of the classical gauge 
generators can provide well-defined quantizations. 

Addressing these questions is crucial, not only for a complete formulation of quantum 
gravity but also for reliable cosmological applications based on the resulting set of equations 
(such as singularity removal or structure formation). In what follows, we will describe the 
current status based on the models of loop quantum cosmology. 



2 Effective dynamics 

In a general sense, effective equations of a quantum system describe the behavior of ex- 
pectation values in a state. Deriving such equations for the expectation values of basic 
operators, such as q and p in quantum mechanics, shows how quantum effects change the 
classical equations of motion. If the equations can be solved or analyzed, the manifestation 
of quantum behavior will be seen. Such effects play an important role in any interacting 
quantum theory, so also in quantum gravity and especially in quantum cosmology which 
devotes itself to the analysis of extremely long evolution times. During those times, quan- 
tum states may change drastically, and quantum corrections grow to be significant. 



2.1 Momentous quantum mechanics 

Effective equations are state dependent since they are equations for expectation values 
in a state with many more independent parameters. In general, quantum fluctuations 
will influence the behavior of expectation values, and must then be included in e f fectiv e 
equations in some form. More generally, and following iBojowald and Skirzewskil (120061 ). 
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we can parameterize a state by all its moments 

q-q p-p 

G a, b = G^^T = {{q _ { g )np _ ^ ))&)weyl (i) 

defined for a + b > 2, in addition to the expectation values (q) and (p). (The subscript 
"Weyl" denotes totally symmetric ordering. We will use the two notations indicated on 
the left interchangeably, at least for small a + b.) For instance, G 2 ' = G qq is the square of 
position fluctuations, and G 1 ' 1 = G qp = \{qp + pq) — (q){p) the covariance. The values of 
moments are not completely arbitrary, most importantly being restricted by uncertainty 
relations such as 

G qq G pp - (G qp f > -j. (2) 

For pure states, the set of moments as defined here is overcomplete; the framework more 
generally allows for mixed states, too. 

All the moments are dynamical. Given a Hamiltonian operator, for every observable O 
we have an equation of motion 

d(6) _ ([6,h]) 

dt ~ ih [6) 

of its expectation value. Specific examples are obtained for the terms in a moment, and 
so we can derive their equations of motion. For the square of position fluctuations, for 
instance, we have 

Introducing Poisson brackets on the space of moments via 

{(A)Am = ^^, (4) 

extended by linearity and the Leibniz rule, equations of motion take Hamiltonian form: 

r\G a ' b 



dt 



{G a ' , H Q } (5) 



with the quantum Hamiltonian Hq := (H). 

In general, the moments are all coupled to one another and to the equations for expec- 
tation values. One can see this by expanding 



H Q ((q),(p),G a > b ) = (H(q,p)) = (H((q) + (q-(q)),{p) + (p-{p)))) 

1 d a+b H((q},( 
a\b\ d(q} a d(p) 



a,b:a+b>2 



(where we assumed the Hamiltonian operator H to be Weyl-ordered in q and p) and 
noticing the coupling terms of expectation values and moments for any non-quadratic 
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potential. The Hamiltonian flow (EJ) or (jSJ) then couples expectation values (q) and (p) to 
all the moments. 

At this stage, we have an exact but usually horribly complicated Hamiltonian descrip- 
tion of quantum evolution. The partial differential equation for a state in Schrodinger's 
formulation is replaced by infinitely many ordinary differential equations for the moments. 
Most of the labor that goes into deriving tractable effective equations consists in extracting 
the required information about expectation values without having to solve for a full state, 
or all its moments, and to specify the regimes where this is reliable. In a semiclassical 
approximation based on near-Gaussian states, for instance, a moment of order a + b is 
typically of the order h^ a+b ^ 2 , giving rise to a natural expansion in powers of fi. To any 
given order, only finitely many moments need be considered^ 



2.2 Harmonic oscillator 

In special systems, equations of motion for the moments decouple automatically to finite 
sets. The best known case is the Harmonic oscillator, whose quantum Hamiltonian is 

H Q = ^ + \m,Htf + \^G« + ±G™ . (7) 

Second-order moments appear, but they do not couple to the expectation values. They 
rather provide the zero-point energy due to quantum fluctuations. Hamiltonian equations 
of motion are only finitely coupled, 

{(p),H Q } = -mu 2 (q) (8) 

{G a '\H Q } = —aG a - l > b+x - mu 2 bG a+hb - 1 . 
m 

For the expectation values we have exactly the classical equations without quantum cor- 
rections, as is well-known for the harmonic oscillator. Here, effective equations for (q) and 
(p) are identical to the classical ones. 

The remaining equations then show how the state evolves once initial values for mo- 
ments have been chosen. For stationary states, for instance, a vanishing covariance G qp 
ensures that fluctuations are time- independent. The covariance is constant in time if 
Qqp _ qvv j m _ muj 2 G qq = 0, or G pp = m 2 uj 2 G qq . With this, the non-classical contribution 
to the energy Hq in (j7|) is mu 2 G qq , where G qq > h/2mu from ([2]). A minimal two-point 

1 Truncating the phase space in this way leads to Poisson manifolds spanned by the moments to a certain 
order. In general, these Poisson structures are degenerate; for instance, there are three independent second 
order moments, forming a space which cannot carry a non-degenerate Poisson structure. Effective equations 
thus make crucial use of Poisson geometry, not sy mplectic geom etry as in other geometric formulations of 
quantum mechanics such as the one going back to lKibbld dl979h . 



m 

dt 
dip) 
dt 

dG a,b 
dt 
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energy \tuxi results for the ground state, derived purely by effective means (even though 
the ground state is not at all semiclassical) . 

Decoupled equations of this form are very useful because they can directly show impor- 
tant aspects of quantum dynamics. Expectation-value equations decouple from the rest 
whenever the algebra of basic operators together with the Hamiltonian is linear: in this 
case, the time derivative of the expectation value of any one of the basic operators is an 
expectation value of a basic operator. Such systems are solvable in a strong sense; there is 
no quantum back-reaction from the moments on the dynamics of expectation values. And 
the dynamics of moments of a given order depends only on moments of the same order. In 
terms of quantum field theory, solvable models correspond to free theories. 



2.3 Low-energy effective potential 

With interactions, or for anharmonic terms, quantum back-reaction results. Moments 
couple non-trivially to expectation values and become important for their dynamics. In 
this way, the state dependence of effective equations ensues. The system of equations 
necessarily becomes of higher dimension than classically, with new dynamical quantum 
degrees of freedom given by the moments. 

For a quantum mechanical system with an anharmonic potential, for instance, we have 
the classical Hamiltonian H = ^p 2 + \moj 2 q 2 + U(q) with anharmonicity U(q). In terms 
of dimensionless quantum variables 

g a,b = h -(a+b)/2( muj y/2-b/2 G a,b (g) 

the quantum Hamiltonian is 

H Q = ^ 2 + lrnu 2 (q) 2 + U m + ^(9 2fi + 9°^ (10) 

n=2 ' ^ ' 



and generates equations of motion 



dt m 
d(p) 



dt 



= —moo 2 



00 1 / h \ n/2 

(?) " U\(q)) " E - — tf (B+1) «?»0 B '° (11) 

n=2 v 7 



g a ' b = awg a - 1 ' b+1 - bug^ 1 *' 1 - bU "^ g ^-i 

muj 



bVhU'"((q}) (^. x _ „ +2)6 _ 1 (6-1) (6- 2) b _ 3 



2{mu) 6 ' z \ 12 
I bhU""((q)) ( _ + (b-l)(b-2) \ + 

6(mui) 2 \ 4 J 

The leading quantum correction appears in the equation for (p) at order h (n — 2), 
for which we have to know the position fluctuation g 2 ' = g qq . In general, this is an 
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independent variable subject to its own equation of motion. Its evolution couples to other 
quantum variables, eventually making the whole infinite system coupled. At this stage, 
approximations are required. For semiclassical states, we may drop terms of higher order 
in providing a closed system of effective equations 

dt m 

^ = -muj\q)-UXm-\ — U'"{(q))9 qq + 0{h^) 
dt I mu 

dg qq 



dt 



2ug 



qp 



dq qp U"((q)) r- 

^- = u{g™-g<* -^Mf + O (VE) 
dt moj 

dt mu 



for expectation values and second-order moments. With higher moments dropped, this 
system shows the leading quantum corrections in a finitely-coupled system. 

For anharmonic oscillators, a further treatment is possible which makes use of an adia- 
batic approximation: we assume that time derivatives of the moments are small compared 
to the other terms. In this case, equations of motion for moments become algebraic rela- 
tionships between them. For the leading adiabatic order g^' 6 , combined with the previous 
^-expansion, we must solve 

= {9 a o\ H Q } = u, (agl^ - b (l + g^' 1 ) + 0(Vh) 

which is of the stationary harmonic form, but with position-dependent coefficients. Al- 
though d(q)/dt ^ 0, this is consistent with the adiabatic assumption. The general solution 
is 



for even a and 6, and g^' = whenever a or b are odd. The values of g^ ' for n even remain 
free, but must satisfy a condition following from the next adiabatic order g\ 



a.b 



For the first adiabatic order, we now consider gQ b weakly time- dependent via (q), but 
assume g a ' b — gQ b time independent. Solutions to the resulting equations {gi' b — go' b , Hq} = 
provide the first adiabatic ordeiEl whose time dependence is given by derivatives of 



2 By definition ([9]), the leading semiclassical //.-dependence of g a,b has been factored out, such that 
explicit factors of h in the equations of motion suffice to read off orders. In terms of the original moments, 
at second order we ignore all terms h n G a,b with 2n + a + b > 2. 

3 The n-th adiabatic order g% b is assumed to satisfy {gn b ~ 9n-i> Hq} = and thus solves {g^' b , Hq} = 
9n-x- Iterating over n, this provides algebraic equations for all orders. 



S 



the zeroth (adiabatic) order moments: 



{g a {\H Q } 



0-1,6+1 

u\ag 1 



Ml 



U"{(q)) \ a+1,6-1 \ _ -a,6 

#1 J — fi'o 



This equation implies 



6 even 



a + b)/2 
b/2 



ma; 2 / 



and requires g^' = C n (l + U"((q})/muj 2 )^ n ^ 4: . The remaining constant C n can be fixed by 
requiring the moments to be of harmonic-oscillator form for U = 0: C n = 2 _n n!/(n/2)!. 
In particular, <?q 9 = 1(1 + U" '((?)) / 'mw 2 )" 1 '' 2 , and the quantum correction to the effective 
force is — \{h/mu)U'''{{q)){l + U" ((q)) / muj 2 )^ 1 / 2 as it arises from an effective potential 



Kff(g) = JmwV + U(q) + ~hu\l+ 1 " (<l) 



(12) 



from ^hu(gQ 9 + g 1 ^) + \{h/muj)U" \(q))g™ in OH]) - This funct ion agrees with path- integral 



calculations of the low-energy effective action as derived by ICametti et all ( 120001 ). Here, 
only the zeroth adi abatic order combined wi t h firs t order in H has been used. To second 
order computed bv iBoiowald and Skirzewskil (1200(1) . there is also a correction for the mass 
term, still in agreement with ICametti et al. teoooh . 

The adiabatic approximation is an example for situations in which the coupled quantum 
evolution of moments can be further reduced to result in explicit effective forces. While the 
forces come from coupling terms between expectation values and quantum variables such as 
fluctuations, the latter do not appear explicitly. Quantum effects then manifest themselves 
indirectly in the form of effective terms depending on the expectation values, their quantum 
origin being indicated only by the presence of h but not by explicit quantum degrees of 
freedom. Such effective terms cannot always be derived since some regimes, where effective 
descriptions may well apply, do require the larger freedom of higher-dimensional effective 
systems coupling quantum degrees of freedom explicitly. For low-energy effective potentials, 
the adiabatic approximation is responsible for the reduction to a system of classical form 
as far as degrees of freedom are concerned. 

From the derivation we can also see how the state dependence is a crucial part of effective 
equations, even though the final expression (112jl for the potential seems state independent. 
To fix all free constants in effective equations for anharmonic systems, in particular C2, we 
had to refer to the harmonic-oscillator ground state. We are thus expanding around the 
vacuum of the free, solvable theory whose moments are known. Quantum corrections from 
the interacting vacuum have been derived in passing, by obtaining the leading adiabatic 
orders of g qq . For the low-energy effective potential, this result was re-inserted in the 
equations of motion, somewhat hiding the state dependence. 

With the higher-dimensional viewpoint of effective systems, keeping some of the mo- 
ments as independent parameters, we obtain extra information about the interacting theory 
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not seen in a simple effective potential. For instance, from our intermediate calculations 
we directly have g q q = §(1 + U" '((<?)) /mw 2 ) 4 / 2 , while g p p = (1 + f/"((g»/mu; 2 )^ 9 = 
|(1 + U" ((q)) / muj 2 ) 1 ! 2 and gg P = 0- To zeroth adiabatic order, the interacting ground 
state keeps saturating the uncertainty relation, but fluctuations are no longer exactly con- 
stant. 



2.4 Quantum cosmology 

In relativistic systems, there is no absolute time with evolution generated by a Hamilto- 
nian. Rather, relativistic systems are subject to a Hamiltonian constraint C. It generates 
arbitrary changes of the time coordinate as gauge transformations S £ f = e{f, C} for phase- 
space functions /. From observable quantities O left unchanged by gauge transformations, 
that is {O, C} = 0, dynamical properties follow. Since the invariance condition {O, C} = 
removes one dimension from the initial phase space, for consistency we must require C = 
as a constraint 

Constrained formulations can be introduced also for non-relativistic systems by pa- 
rameterization, adding a time degree of freedom t with momentum p t and replacing the 
Hamiltonian H by the Hamiltonian constraint C = p t — H. For time-independent Hamil- 
tonians, p t is gauge invariant while the gauge transformation of t is S e t = e{t, C} = e; t 
can be changed at will. For the remaining observables, 

= {O 1 C} = ^-{O,H} 

imposes Hamilton's equations of motion. 

For relativistic quantum systems, the effective techniques described so far cannot di- 
rectly be app lied. There are extensi o ns to effective constraints briefly described later, as 



developed by iBojowald et all (l2009al ) ; iBojowald and Tsobanjanl (120091 ) . But more simply 



effective equation techniques can be used if systems are first deparameterized, reverting the 
above procedure. If there is a variable 0, then called internal time, such that the Hamil- 
tonian constraint can be written as C = p^ — H(q,p) with the momentum p^ of and a 
function H independent of and p^, gauge transformations generated by the constraint 
take the form 

Sef = e{f(q, p, 0, p,), C} = e(^-- {/, H} 



4 Factoring out the Hamiltonian flow generated by the constraint C via its Hamiltonian vector field 
Xc = {'jC}, we obtain a projection it: M — > M/Xc from the original phase space M by identifying all 
points along the orbits of Xc- All observables O naturally descend to the factor space since they are 
constant along the orbits, and so does C. In this way, we obtain a complete set of functions on M/Xc- 
On the factor space, we have a natural Poisson structure {/,s}m/x c = I 71 "*/' 7r *.9}-A/i pushing forward 
the Poisson bivector via tt. This Poisson structure is degenerate: {O, C}m/x c = f° r au functions O 
on M/Xc- The constraint C becomes a Casimir function on the factor space, and symplectic leaves of 
the Poisson structure are given by C = const. Any leaf carries a non-degenerate symplectic structure 
and can be taken as a reduced phase space, but C = is distinguished: In this case, we can write gauge 
transformations as 6 e f — {/, eC} even for phase-sp ace functions e. More detai ls of Poisson geometry in 



the context of constrained systems are described bv lBoiowald and Strobll (J2003). 
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Gauge invariant quantities of the theory are thus those evolving in the usual Hamiltonian 
way as generated by the 0-Hamiltonian H . 

Once deparameterized, the observables of a constrained system can be derived by ana- 
lyzing an ordinary Hamiltonian flow. At this stage, effective techniques as described before 
can be applied to quantizations of deparameterized models. Effective equations of motion 
derived from 

d<6> ([d,H]) 



ih 



(13) 



then provide means to solve for quantum observables (O) (0) and their physical evolution. 
For the initial constrained system, solutions (O)(0) are observables as functions on the 
full phase space including (0,p</,). (Remaining phase space variables enter the expression 
(O)(0) via initial values taken for them when solving the differential equations ( fl3|) .) The 
deperameterization endows (0)((p) with the interpretation as an observable (O) evolv- 
ing with respect to the internal time <fi. This is the rela t ional p icture f o r inter preting 
constrained dynamic s, developed classically by iBergmann (119611 ). iRovellil (ll991bU ah and 
Dittrichl f l2007l . l2006h . 

Examples for such systems in cosmology are homogeneous models sourced by a free, 
massless scalar. Its energy density is purely kinetic, p = ^a~ 6 p^, such that the Fried- 
mann equation, solved for p^, provides a Hamiltonian for ^-evolution. For a spatially flat 
Friedmann-Robertson-Walker model, the 0-Hamiltonian then turns out to be quadratic 
in suitable canonical variables: for instance, the Hubble parameter T-L = a /a is canonically 
conjugate to the volume V = a 3 , and the Friedmann equation 



8ttG 



4ttG: 



tells us that p 2 , oc V 2 l-i 2 . Upon taking a square root, we have a quad ratic ^-Hamiltoni an. 



Based on this observation, state properties have been determined by iBojowaldl (j2007bU al). 
(Strictly speaking, the 0-Hamiltonian is of the form \qp\ which is not quadratic. However, 
for effective equations one can show that the absolute value can be dropped, providing 
a linear quantum system. We only have to require an initial state to be supported on a 
definite part of the spectrum of qp, either the positive or the negative one, which is then 
preserved in time since H is preserved. The absolute value just amounts to multiplication 
with ±1. For initially semiclassical states, this requirement does not lead to restrictions 
for low-order moments.) 

For other systems, perturbation theory can be used as described above f or anh armonic 
oscillators. A crucial difference, according to the analysis by lBojowald et al.l (j2007al ). is that 
no adiabatic regime has so far been found for quantum cosmology, blocking the complete 
expression of quantum variables in terms of an effective potential. On the other hand, 
higher- dimensional effective systems, where quantum variables are taken as independent 
variables subject to their own evolution, can be analyzed and show how states back-react 
on the expectation-value trajectories. 
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If we include a mass term or a potential for the scalar, the system becomes "time- 
dependent" in 0. Extra care is required, but perturbation theory still applies for small 
and flat potentials. The justification for this procedure in the time-dependent case comes 
fro m an extension of the e f fective-equation procedure to co nstrained systems introduced 
by Bojowald et al.l ( 2009a ): Bojowald and Tsobanjan ( 2009 ) without requiring deparam- 
eterization. (Quantum constraint operators then imply the existence of infinitely many 
constraints (polC) on the quantum phase space, in general all independent for different 
polynomials pol in the basic operators. This large number of constraints restricts not only 
expectation values (as the classical analog) but also the corresponding quantum variables. 
A complete reduction to the physical state space results, without requiring a deparameter- 
ization to exist. With these general techniques, effective equations for relativistic systems 
are thus fully justified. As discussed in Sec. HI inhomogeneous models bring in a new level 
due to the anomaly problem, which does not arise for a single classical constraint. 



2.5 Symplectic structure 

Quantum corrections in canonical effective equations come from corrected Hamiltonians 
or corrected constraints as expectation values of Hamiltonian (constraint) operators. The 
Hamiltonian H or constraints Cj, together with Lagrange multipliers N 1 (for gravity lapse 
and shift multiplying the Hamiltonian and diffeomorphism constraints), form an important 
contribution to the action. But the action 

S[q(t),p(t);N I } = J dt (qp-Hiqrf-rfdfap)) 

of a Hamiltonian system in canonical form has an extra contribution, the one that de- 
termines the symplectic structure between configuration and momentum variables. One 
might wonder whether it is enough to look for quantum corrections in the Hamiltonian 
or the constraints without correcting the symplectic structure. With symplectic-structure 
corrections, an effective action might take a different form than suggested by an analysis 
only of Hamiltonians and constraints. 

There might indeed be corrections to the symplectic structure, but they would follow 
from the same algebraic notions used for canonical effective equations. Poisson brackets 
of quantum variables are given by expectation values of commutators as used before; see 
OH). Any potential corrections to the symplectic structure have thus been taken care of 
by applying that formula consistently. If there are no changes in Poisson-bracket relations 
for effective equations, it is only because they do not change for expectation values of 
basic operators: they have commutator relations mimicking the classical Poisson algebra, 
which is linear for the basic objects. Taking expectation values of these linear structures 
does not lead to corrections. For canonical basic operators q and p, for instance, we have 

{<g>,(p)} = ([g,p])M=l. 

The Poisson structure changes only because the dimension of the phase space increases 
by the quantum variables G a,b . These variables satisfy Poisson relationships following from 
the quantum theory, but this does not affect the Poisson brackets of expectation values of 
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basic operators. In an action, the symplectic term for classical variables remains unchanged. 
But if a higher-dimensional effective system with independent quantum variables is used, 
their symplectic terms add to the action. The symplectic structure is only extended to 
include new degrees of freedom; it is not quantum corrected. 

Sometimes, one can make assumptions^ about the dependence of quantum variables on 
expectation values, or even derive those by an adiabatic approximation. If this is done, one 
can insert such expressions, schematically G I (q,p), into the symplectic form for quantum 
variables, of the form Qjj(G)dG I AdG J where Qjj(G) follows from (j4j) purely kinematically. 
Then, a term 2Qu{dG I /dq)(dG J /dp)dq A dp results which would add to the symplectic 
term of expectation values. Now, the correction has dynamical information via the partial 
solutions G (q,p). (For the specific example of anharmonic oscillators, no such corrections 
arise since the effective equations for g ab , and thus their solutions, depend only on (q), not 
(p).) With independent quantum variables, as they are most often required, no symplectic- 
structure corrections result. Expectation values of Hamiltonians or constraints are then 
the key source for quantum corrections. 

For quantum gravity, we expect higher-curvature terms in an effective action, and thus 
higher- derivative terms. By the preceding discussion, this can only come from independent 
moments G a,b . But quantum gravity has additional implications, such as the emergence of 
discrete spatial structures expected in particular from loop quantum gravity. They, too, 
must affect the effective dynamics. 

3 Discrete dynamics 

Our main interest from now on will be to apply effective equation techniques to canonical 
quantum gravity. A strict derivation requires detailed knowledge of the mathematical 
properties of operators involved, as well as information about the representation of states. 
Effective equations and effective constraints are, after all, obtained via expectation values 
of Hamilonians or Hamiltonian constraint operators. These operators are first constructed 
within the full quantum theory of gravity or one of its models, and the resulting constructs 
will be sensitive to properties of the basic operators, analogous to q and p in quantum 
mechanics. Some of these properties will descend to the effective level once expectation 
values of Hamiltonians are computed, and then affect also the effective dynamics. At this 
stage, general effective techniques naturally tie in with specific constructions of a concrete 
quantum theory at hand. 

3.1 Loop quantum gravity 

Classically, the canonical structure of general relativity is given by the spatial part q a h 
of the space-time metric g a b, as well as momenta related to its change in time, or the 
extrinsic curvature K a i, of the spatial slices E. These specific quantities are meaningful 

5 A quantum cosmological model with assumptions for semiclassical states has been analyzed bv lTaveras 
(|2008h . 
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only when a choice for time, a time function t such that E: t = const, has been made, and 
so this formalism is often seen as breaking covariance. But space-time covariance is broken 
only superficially and is restored when all the dynamical constraints have been solved - 
in addition to a local Hamiltonian constraint, three components of the diffeomorphism 
constraint to generate all four independent space-time coordinate changes. The theory 
is, after all, equivalent to general relativity in its Lagrangian form; just formulating it 
in different variables cannot destroy underlying symmetries. In fact, the Hamiltonian 
formulation still shows the full generators of all gauge transformations in explicit form, 
by the constraints it implies for the fields. By the general theory of constraint analysis, 
a discussion of gauge at the Hamiltonian levels is then even more powerful than at the 
Lagrangian one. 

3.1.1 Smearing 

Before addressing quantum constraints, the quantum theory must be set up. For a well- 
defined quantization one turns the basic fields into operators such that the classical Pois- 
son bracket is reflected in commutator relationships. In constructions of quantum field 
theories one has to face the problem that Poisson brackets of the fields involve delta 
functions since they are non- vanishing only when the values of two conjugate fields are 
taken at the same point, as in {(j)(x) , p</,(y)} = S(x,y) for a scalar field. A simple but 
powerful remedy is to "smear" the fields by integrating them against test functions over 
space. Again for a scalar field, we could use (f>\p] := f d 3 xy / dlit^ / u(:r)0(a;) for which 
{4>[lAiP(j>{y)} = f d 3 xfi(x)5(x,y) = fi(y). The Poisson algebra for the smeared fields is 
free of delta functions thanks to the integration. Moreover, for sufficiently general classes 
of test functions n(x), smeared fields capture the full information contained in the local 
fields and can be used to set up a general theory. 

After smearing, a well-defined Poisson algebra of basic objects results, ready to be 
turned into an operator algebra by investigating its representations. For gravity, we would 
use smeared versions of the spatial metric and its change in time. But here, a second 
problem arises. One of the dynamical fields to be quantized is the spatial metric, but to 
smear fields we need a metric for the integration measure. This is no problem when fields 
to be quantized are non-gravitational. As with the scalar in the example above, we would 
use the background space-time on which the scalar moves, obtaining quantum field theory 
on a given, possibly curved background. But what do we do if the metric itself is one of 
the fields to be quantized? If we use it to smear itself, the resulting object becomes ugly, 
non-linear and too complicated as one of the basic quantities of a quantum field theory. 
If we introduce a separate metric just for the purpose of smearing the dynamical fields 
of gravity, this auxiliary input is likely to remain in the results derived from the theory. 
We would be quantizing gravitational excitations on a given space-time background, not 
the full gravitational field or space-time itself. We would be violating the great insight of 
general relativity by formulating physics on an auxiliary space-time, rather than realizing 

6 For an introduction to these techniques see iBoiowaldl (|2010f ). 
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gravity as the manifestation of space-time geometry. 

Fortunately, it is possible to smear fields and yet avoid the introduction of auxiliary 
metrics. Back to the scalar, we can choose to smear p^ instead of 0. The momentum 
of a scalar field is a scalar density; it transforms under coordinate changes with an extra 
factor of the Jacobian for the coordinate transformation. This is a direct consequence 
of the definition p^ = dC/dcf) as a derivative of the Lagrangian density. Explicitly, the 
canonical variable p^ = y/ det q<fi already carries the correct measure factor which need not 
be introduced by an auxiliary metricJ3 The smeared version p^[A] := / d 3 yX(y)p^,(y) is 
well-defined for any function A, and it suffices to remove delta functions from the Poisson 
algebra: {<p(x) , p^[X]} = X(x). 



3.1.2 Holonomies and fluxes 



For tensorial fields as we have them in gravity, background-independent smearings are often 
more difficult to find. Loop quantum gravity has provided suitable procedures for general 
relativity, but for this it must first transform from metric variables to connections with 
their conjugates, densitized triads. Connections and densitized vector fields turn out to 
have just the right transformation properties under coordinate changes that they can, with 
one loopy trick, be smeared background independently. A well-defined quantization results 
with several immediate implications for the basic operators encoding spatial geometry, as 
well as far-reaching and sometimes surprising consequences in the resulting dynamics. 

Instead of using the spatial metric q a b, spatial geometry is expressed by a densitized triad 
Ef = y/d^qe^ such that EfE\ = det qq ab . The densitized triad is canonically conjugate 
to K % „ := Kab^ 1 in terms of extrinsic curvature K a h. T o obtain a connection with its 



useful transformation properties, we finally follow lAshtekarl (119871 ) and iBarbero G.I ( 119951 ) 
to introduce the Ashtekar-Barbero connection A\ = T l a + with the spin connection T l a 
compatible with the densitized triad and a positive real num ber 7, the Bar bero-Immirzi 
parameter (whose role for quantum geometry was realized by llmmirzil (119971 ) ) . 

The elementary objects loop quantum gravity takes for its representation are holonomies 
and fluxes, 



h e (A) =Vexp( [ dse a AirA and F S (E) = [ d 2 yn a Efr t , 



first used in this context by iRovelli and Smolinl (ll990l ). A key advantage is that their 
algebra under Poisson brackets is well-defined, free of delta-functions (unlike the algebra 
of fields), and yet independent of any background metric. Only the dynamical fiels A\ and 
Ef are used, together with kinematical objects such as curves e C E and surfaces ScE 
as well as the tangent vectors e a and co-normals n a to them, but no independent metric 
structured All spatial geometrical properties are reconstructed from Ef via fluxes, and 



7 Although y'det q appears in the relationship between p^ and </>, from the viewpoint of Poisson geometry 
p^ (but not (/>) is independent of the metric: {p<j>,p ab } = for the momenta p ab of q a b- 

8 The co-normal, unlike the normal n a , is metric independent: for a surface S: f = const, n a — (df) a - 
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space-time geometry follows with A l a via holonomies once equations of motion (or rather 
the constraints of relativity) are imposed. In this way, loop quantum gravity provides a 
framework for background-independent quantum theories of gravity. 

Once a well-defined algebra of basic objects has been chosen, one can determine its 
representations to arrive at possible quantum theories. In the connection representation, a 
complete set of states i/){A % a ) is generated by holonomies as multiplication operators acting 
on ip(A l a ) = 1. In the case of loop quantum gravity, this has an immediate and general 
consequence. Holonomies take values in SU(2), and fluxes, depending on the momenta 
Ef, become derivative operators on SU(2) — just like angular momentum in quantum 
mechanics. For a dense set of states, only a finite number of holonomies (along curves 
intersecting the flux surface S) contribute to a given flux. Finite sums of angular mo- 
mentum operators with discrete spectra have a discrete spectrum, too: spatial geometry 
is discrete; flux operators quantizing the densitized triad and thus encoding spatial ge- 
ometry acquire discrete sp ectra. So do spatia l geom e trical quantities such as are a s and 
volum es as constructed by iRovelli and Smolinl (119951 ); lAshtekar and Lewandowskil (11997 . 
19981 ). No extra assumptions are required; one merely has to fix the basic algebra and 
follow mathematical procedures to analyze its representations. A different algebra might 
lead to other properties, possibly not with discrete spatial geometry. But no alternative 
procedure providing a well-defined and smeared, yet background independent quantization 
has been found. The holonomy-flux algebra suggested by its natural smearing, on the 
other han d, has a unique irreducibl e, c yclic, diffeomorphis m covariant representation as 
proven by Lewandowski et al. ( 2006 ) an dFleischhack ( 2009 ). Most of these properties are 
described by ISahlmannl ( 120101 ) . 



3.1.3 Dynamics 

Kinematical properties are elegant, simple, and largely unique. The theory starts to get 
considerably more messy when its dynamics is considered. Here, two major tasks must 
be performed: Dynamical operators, mainly the Hamiltonian constraint, must be defined 
from the basic ones, holonomies and fluxes. This is the constructive part of the task. For 
suitable constructions of the Hamiltonian constraint, the dynamics of the theory must then 
be evaluated, a process still full of several open issues. 

At the constructive stage, many choices can be made, and strong consistency condi- 
tions must be respected. We are witnessing an epic battle between the liberating anarchy of 
choice and the uniformizing tyranny of constraints. Just writing a Hamiltonian constraint 
operator is tedious but possible in many ways. There are ubiquitous factor ordering am- 
biguities, as well as other choices specific to loop quantum gravity. On the other hand, 
the Hamiltonian constraint provides not only an equation of motion, it also generates a 
crucial part of the gauge transformations responsible for general covariance. A consistent 
quantization must keep gauge degrees of freedom as gauge, and not overly restrict the 
number of physical, non-gauge degrees of freedom. All this can usefully be formalized, as 
we will see later, providing strong algebraic conditions. They are so strong that to date, 
despite the abundant freedom of choices in constructing Hamiltonian constraint operators, 
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no full consistent version has been found. 

Once a consistent version of the dynamics exists, it must be evaluated. We must 
find solutions, and determine the physical observables they provide. Their values, finally, 
can be used for predictions such as small deviations from the expected classical behavior. 
Since no full consistent version has been found yet, and since all potential condidates are 
highly complicated, construction issues of the dynamics have so far dominated strongly over 
evaluation issues. In several model systems, on the other hand, the dynamics can often be 
simplified so much that it can be analyzed rather explicitly. Many useful techniques are 
now available, mainly in the context of effective equations. We will come back to these 
applications after discussing more details of the construction side of the problem. 



3.1.4 Hamiltonian constraint 

Specifically, the Hamiltonian constraint of general relativity in Ashtekar variables is 



C[N) = Jd 3 xN(x) U jk Fl b - 2(1 + i- 2 )KKi^ E ^- 



v/|det£| Vl det ^l 

with the curvature F l ah of the Ashtekar connection and extrinsic curvature K % a . It has to 
vanish for all lapse functions N(x), thus providing infinitely many constraints. If C[N] 
is to be turned into an operator, using the basic expressions for holonomies and fluxes, 
several obstacles must be overcome. 

First, there is the potentially singular inverse determinant of the densitized triad. No di- 
rect quantization exists since the densitized triad has been quantized to flux operators with 
discrete spectra, containing zero. Such operators lack densely defined inverses. Neverthe- 
less, quantizations with the appropriate inverse as the semiclassical limit can be obtained 
making use of the classical identity 

{^/vfdSiid-,}*^.-^ (") 



or variants) as introduced by iThiemannl (11998al ). There is no inverse on the left-hand 



side. Instead, the expression involving the densitized triad is the spatial volume which 
(with some regularization) can directly be quantized. The connection components can be 
expressed via holonomies, and the Poisson bracket will, at the quantum level, be quantized 
to a commutator divided by ih. A well-defined operator results with the right-hand side 
of ( jT4"j) as the desired semiclassical limit by construction. 

The remaining factors in the Hamiltonian constraint involve the Ashtekar curvature as 
well as extrinsic curvature. For the curvature components F* b , we can use 

4^(^ = A^(^a-1) + 0(A) (15) 

where A is a small loop starting at a point x, spanning a coordinate area A, and with 
tangent vectors s" and at x. On the right-hand side, the holonomy h\ can readily 
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Figure 1: Schematically, the local action of the Hamiltonian constraint operator on a 
state. States are generated by holonomies as multiplication operators, visualized by the 
graph formed by all curves e used. Moreover, labels j e determine matrix representations 
of SU(2)-valued h e (A). Due to ( 1T5|) . new curves and vertices are typically created when 
the Hamiltonian constraint acts. 



be quantized, and to leading order provides curvature components as required for the 
constraint. 

Extrinsic curvature, finally, is a more complicated object in terms of the basic ones but 
can be obtained from what has been provided so far: 



, , e*> k EfE k 



b 



K a oc < A a , < / d 6 xF l ab , J =, / x/| detEWx 



^/\detE\ , 



expresses extrinsic curvature in terms of a nested Poisson bracket involving the spatial 
volume and the first term of the constraint, already provided by the preceding steps. 

In this way, holonomy and flux operators make up the Hamiltonian constraint oper- 
ator C as a densely defined operator (includ ing, in the non- va cuum case, regular matter 



Hamiltonians again exploiting (I14p following iThiemannl ( I1998bl )). It determines the phys- 



ical solution space by its kernel: physical states if} (A), assumed again in the connection 
representation, must satisfy Cif>(A) = 0, or (C + 87rGi^ matter )^(A, 0, . . .) = if matter is 
present. The action of the constraint is rather complicated, as visualized schematically in 

Fig. m 

Physical states of interest normally belong to zero in the continuous part of the spectrum 
of C, which requires the introduction of a new physical Hilbert space spanned by the 
solutions to the quantum constraint equation and equipped with a suitable physical inner 
product. Constructing the physical Hilbert space leads over to evaluating the dynamics, 
for physical states would provide predictions from expectation values of observables. But 
this stage has been brought to completion only in a few simple (and very special) models; 
in general the outlook toward a full implementation is rather pessimistic. 

Here, effective constraint techniques become useful because they allow one to address 
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physical properties, corresponding to observables in the physical Hilbert space without hav- 
ing to deal explicitly with states. Conditions from the physical inner product are rather 
implemented by reality conditions for expectation values and quantum variables such as 
fluctuations, which can be done more simply and more generally than for entire states. Ef- 
fective constraints thus provide a good handle on generic properties of physical observables, 
at least in semiclassical regimes where not all the moments need be considered. Especially 
in cosmological situations they provide an ideal framework. In inhomogeneous contexts, 
they allow a detailed discussion of the anomaly issue, and show whether the different effects 
expected from the basic operators of quantum gravity can lead to a consistent form of the 
dynamics. 

Consistency conditions in the presence of ambiguities are useful, for they constrain the 
choices. But the question remains whether loop quantum gravity can be fully consistent 
at all — implementing covariance at the quantum level might, after all, not leave any 
consistent physical states. If the consistency conditions are weak, on the other hand, 
ambiguities would remain even at the physical level. In between, at a razor-thin balance 
between ambiguity and constraints, lies the case of a unique consistent theory, a possibility 
which exists in loop quantum gravity but for which at present no evidence has been found. 

In model systems it has at least been shown that consistency in the form of covariance 
can be achieved even in the presence of quantum corrections resulting from the discreteness. 
This statement may come as a surprise, for discrete spatial or space-time structures are 
naively expected to break local Lorentz or even rotational symmetries. Nevertheless, at 
the level of effective equations one can see that covariance can be respected - - but it 
cannot leave the classical algebra of local symmetries invariant. While quantum-gravity 
corrections provide a consistent deformation of general relativity — preserving the number 
of independent gauge generators — the algebra is truly deformed. Gauge transformations 
are no longer local Lorentz or coordinate transformations, but of a different type relevant 
for the quantum space-time structures derived from loop quantum gravity. 

Much of these conclusions makes use of results obtained through several years for models 
of loop quantum cosmology. We will now start a general exposition beginning with the 
simplest isotropic models, introducing the key features (and sources of ambiguities), and 
then leading over to the discussion of covariance in the next section. 



3.2 Isotropic loop quantum cosmology 



Loop quantum cosmology provides quantizations of symmetry-reduced models of general 
relativity, starting with isotropic ones in the simplest cases. It is thus a minisuperspace 
quantization, but not only that. There are explicit relationships between the states and the 
algebra of basic operators in a model of loop quantum cosmology, which shows how they 
descend from the analogous expressions in the full theory. For different kinds of reductions 
such as im plementing; isotropy, h o mogeneity or spherica l symmetry, this has been con- 
structed by lBojowald et al.l ( 120061 ); iBojowaldl ( 12004 120061 ). In particular, properties of the 
basic operators, for instance the discreteness of the flux spectra, are preserved and realized 
also in cosmological models. Crucial implications for the discreteness of spatial geometry 
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can then be analyzed in their dynamical context. Qualitatively, all applications of loop 
quantum cosmology rel y on this preservation of discre teness properties by the symmetry 



reduction introduced by iBojowald and Kastrupl (120001 ) 



While existing methods do not allow one to derive the dynamics of models directly from 
full constraint operators, it can be constructed in an analogous way using all the construc- 
tion steps sketched for the full constraint. Also here, crucial properties are preserved. And 
even though the incompletely known relationship to the full theory introduces additional 
ambiguities in the dynamics of models, reliable conclusions can still be drawn provided the 
considerations are generic enough. Here, the need for sufficiently general parameterizations 
of ambiguities arises. 

3.2.1 Basic variables 

In spatially flat isotropic models, the basic canonical fields reduce to A\ = c8 l a and Ef = P5" 1 
with two dynamical variables c and p. Relating the densitized triad to the spatial metric 
shows that \p\ = a 2 is given by the scale factor. (Due to the freedom of choosing an 
orientation of the triad, p can take both signs.) Classically, c = 7a for spatially flat 
models. 

To define the Poisson structure, we pull back the full symplectic form (SivyG)^ 1 f d 3 x5A l a A 
5E°- by the embedding (c,p) (A l a ,E°-) = (cS^pSf). Due to homogeneity, the resulting 
integral diverges if we integrate over all of space and space is infinite, but homogeneity 
also implies that we may just integrate over any finite chunk of coordinate volume Vo, say, 
and still get the complete reduced symplectic form, 3Vo(8ii'yG)~ 1 dc A dp. It depends on 
the arbitrary Vo, which can be hidden in the canonical variables by redefining them as 

c = Vq^ 3 c and p = V^ 3 p . 

In the end, physical quantities must be ensured to depend only on combinations of c, p 
and possibly other ingredients such that they are insensitive to changes of Vq. (Changing 
Vo is independent of changing coordinates, thereby rescaling the scale factor. While c and 
p depend on the scaling but not on Vo, c and p depend on V but not on the scaling.) 

Isotropic connections and densitized triads result in specific versions of holonomies 
and fluxes. They take especially simple forms when evaluated along curves or surfaces 
making use of structures provided by the homogeneity group, i.e. curves along generators of 
translations and surfaces transversal to them. In a homogeneous context, fluxes are simply 
the triad components multiplied with the coordinate area of the surface, and holonomies 

1 /3 

can be computed explicitly: h ej (A l a ) = exp(/icr 7 ) with p = io/V for a curve tj of 
coordinate length £q along the direction e°j = {d/dx^) a in Cartesian coordinates. If we 
take surfaces for fluxes of edge length the same size as ej, fluxes are Fs^E" 1 ) = p 2 pTj for a 
surface Sj transversal to {d/dx^) a . 

These curves and surfaces are particularly useful because they can be arranged in a 
regular lattice of spacing £ . While such a setting is not the most general one (and other 
curves would lead to differen t expr essions for holonomies as computed for instance by 



Brunnemann and Fleischhackl (120071 )). it allows one to capture all the degrees of freedom 
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in a homogeneous model. More importantly, it allows for rough contact with the full 
theory, for the regular lattices envisioned here can be thought of as lattices for states 
in the full theory In particular for comparisons of actions of the dynamical operators, 
such a relationship is convenient. It also allows us to find a useful interpretation of the 
parameter [i characterizing isotropic holonomies and fluxes: For a regular lattice of spacing 
to, contained in a region of size V , /i 3 = il/V =: Af^ 1 is the inverse number of lattice 
sites. This parameter can be seen as providing information about an underlying discrete 
inhomogeneous state which appears even in a homogeneous reduction. In fact, properties 
of M such as its size or possible dynamical features play a crucial role in Hamiltonian 
quantum evolution. 



3.2.2 Representation 

From isotropic holonomies and fluxes, we can construct their loop representation by analogy 
with the full theory. States in the connection representation are functionals of holonomies 
and are thus superpositions of the basic states or exp(i/ic) as functions of the isotropic 
connection component. Here, /i is an arbitary real parameter. An inner product of states 
ca n be derived from integ ration theory on spaces of connections, as developed for instance 



by lAshtekar et al.l ( 119951 ). For isotropic models, this makes all states orthonormal. 

Upon completion to a Hilbert space, a general state takes the form of a countable 
superposition of exp(z/ic) with fi e R: 

^( c ) = Yl /jexp(z^jc) (16) 

IeXCR, countable 

such that \fi\ 2 exists- These states form a Hilbert space equivalent to the £ 2 -space 

formed by the normalizable sequences (fi)i<zx fo r all countable I C K. In this form, 
isotropic states directly result from the reduction of a full state, where all the holonomies 
reduce to exponentials of c with different exponents. Alternatively, the Hilbert space may 
be characterized as the space of square integrable functions ip(c) on the Bohr compactifi- 
cation of the real line, a compact space containing the real line densely and equipped with 
the measure 

dp(c)ip(c) = lim — / ^(c)dc 

C^oo J_ c 

using the normal Lebesgue measure on the right-hand side. This Hilbert space is non- 
separable. 

Basic operators act by multiplication or differentiation: 

exp(i<5c)|/i) = |/i + 5) (17) 

m = y74/^> (18) 

as it indeed follows from the unique holonomy-flux algebra in the full theory as an induced 
representation. A Wheeler-DeWitt representation, by contrast, would not be related to 
such a representation of the full theory. 
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Properties of the loop representation are markedly different from the Wheeler-DeWitt 
one, but they closely mimick properties of the full holonomy-flux representation: 



There is a discrete spectrum of p. While there is a continuous range for /x, all 
eigenstates \fi) are normalizable. For a non-separable Hilbert space as we are dealing 
with here, normalizability of eigenstates does not imply that the eigenvalues form a 
countable subset of the real line. In such a situation, the normalizability condition is 
more general, for it is insensitive to what topology one uses on the set of eigenvalues. 
Even if any real number can appear as an eigenvalue, this set would be discrete if 
one uses a discrete topology of the real line (for instance, one where every subset is 
an open neighborhood). Normalizability of eigenstates will also be one of the crucial 
properties for consequences of flux spectra in loop quantum cosmology. 

There is no operator for c, and only holonomies are represented. Trying to derive an 
operator for c from holonomies, for instance by taking a derivative of the action of 
exp(i5c) by 5 at 5 = 0, fails because holonomies are not represented continuously in 
5: (ju|exp(i5c)|ju) = (fi\fi + S) = 5s,o is not continuous. 



Ashtekar et al.l (120031 ) and lFewster and Sahlmannl (120081 ) discuss these and other properties 



of the Bohr compactification as used in loop quantum cosmology, and IHusain and Winkler 



(12004J ) do so in a quantum cosmology based on ADM variables. 

Since these are the same basic properties as they are realized for the full holonomy- 
flux algebra, they have the same qualitative implications for the dynamics. Once used 
for the construction of operators such as the Hamiltonian constraint, they lead to specific 
quantum-geometry corrections in loop quantum cosmology as they do in the full theory. 
Specifically, 

1. Only almost-periodic functions of c are represented as operators on states (Tl6|) . and 
must be expressible as ip(c) = J2ieicR countable fi exp(zJjc). The Hilbert space does 
not allow an action of c on its dense subset of triad eigenstates; any appearance of 
c not of almost-periodic form, such as the polynomial in the isotropic Hamiltonian 
constraint, must be expressed in terms of almost-periodic functions by adding suit- 
able higher-order corrections in c. They become significant for large values of the 
curvature, via Iqc = c/J\f 1 ^ 3 for a regular distribution of edges where 5i ~ A/" -1 / 3 . 

2. The isotropic flux operator p has a discrete spectrum containing zero and so lacks a 
direct, densely defined inverse. Well-defined versions can be obtained via identities 
such as 

l -e^{e- iSc , \p\^} = {c, \p\^} = ^b^- 1 sgnp (19) 



which mimick the crucial one (JT4J) used in the full theory. For < r < 2 we are 
expressing an inverse of p on the right-hand side, but do not need an inverse on 
the left-hand side; well -defined a n d even bounded operators for inverse powers of 



p result, as derived by iBojowaldl (j2001bl ). Within this range, r is unrestricted by 
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general considerations; it thus appears as an ambiguity parameter (see iBojowald 



(l2002bl ) for further discussions of ambiguities). Also here, corrections to classical 
expressions arise, in this case for small flux values l\p = p/M 2 ^ 3 near the Planck 
scale. From a quantization of the left-hand side of (Tl9l) . eigenvalues can readily be 
derived. They have the form |7T7Gr|p| r / 2 ~ 1 a r (p) where a r (p) ~ 1 asymptotically for 
large p/N" 2 ^ 3 3> £p while a r (p) — » rapidly for p — > 0, cutting off the divergence of 
|p|r/2-i_ rj ne correc ^ semiclassical limit is guaranteed by a r approaching one, which 
it does from above irrespective of quantization ambiguities. 

3.2.3 Difference equation 

The main quantity for which these considerations play a role is the Hamiltonian constraint, 
in isotropic variables providing the Friedmann equation 



c 2 ^\p\ SttG 

= — - — n. 
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2 



matter 



with the matter Hamiltonian i? ma tter — Voa 3 p. Classically, its dependence on c is via c 2 , 
which is not almost periodic. There are many ways to express c 2 in terms of almost- 
periodic functions such that they approximate c 2 for small curvature, c < 1. In general 
terms, properties of the quantum representation space require the replacement of c 2 to be a 
bounded function of c, given by a normalizable superposition exp(5/c) with 5j in a countable 
subset of the real line. Each exponenti al acts by a shift of triad eigenvalues, rather than 
a derivative operator, and so, following iBoiowaldl fl200ld . l2002at ). upon loop quantization 



the Hamiltonian constraint equation becomes a difference equation for components of a 
wave function in the triad representation. 

A choice usually made is c 2 ~ 5~ 2 sin 2 (5c), but others are possible, constituting further 
ambiguities. A Hamiltonian constraint operator 

3 - — - 2 ^7 ==z - 

Ciso = Sinfc VM - -^matter 

87rG7^ 

results. Expanding a state = J^VvIa 4 ) in the triad eigenbasis Asol^) = is 

equivalent to the difference equation 

CO + 25)^+ 25 - 2CO)W + CO - 25)^-25 = ^7 2 <^matter(/^ (20) 



as derived by iBojowaldl (l2002al ). where C(n) are eigenvalues of y/\p\, e.g. proportional to 
\f\p\- Additional matter fields have been suppressed in the notation, which would provide 
further independent variables in ip^ acted on by the matter Hamiltonian -ff ma tter(/-0 • (Unless 
there are curvature couplings, matter Hamiltonians depend on the densitized triad but not 
on the connection; the right-hand side of (1201) is then not a difference expression. Fermions, 
which require a coupling to the gravitational connection, have been discussed in this context 
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bvlBoiowald and Dad (120081 ) , and non-minimally coupled scalars by IBojowald and Kagan 
d2006h .) 

Due to the inequivalence of representations, this difference equation replaces the differ- 
ential Wheeler-DeWitt equation. At this dynamical level, after choosing the representation 
of c 2 by almost-periodic functions, the dynamics can be restricted to a separable subsector 
of the initial non-separable Hilbert space. From this perspective, one could have started 
directly from a separable Hilbert space by choosing a specific countable set {hi} C R, such 
as Hi = Ifi with integer / as originally defined by IBojowald! (j2002al ) . Many dynamical 
properties would follow in the same way since non-separable features appear only at the 
kinematical level. 

However, there is a good reason for dealing with a non-separable kinematical Hilbert 
space, allowing for all real values of /i: To capture the most general and viable dynamics, 
one expects lattice refinement to happen, with a discreteness scale changing dynamically. 
Indeed, full Hamiltonian constraint operators create new edges and vertices, changing the 
graph on which a state is defined, as illustrated in Fig. [TJ In an isotropic context, a trace of 
this feature must still be left since the dynamical parameter 5 came from the edge length 

1 /3 

used in elementary holonomies. While Vo in 5 — £o/V remains constant as a mere auxil- 
iary parameter introduced in the setup, £o must be adapted to the geometry. Geometrical 
distances, after all, increase as the universe expands and the region of coordinate size Vo 
has a growing geometrical size Voa 3 . If a lattice state for a larger region has more sites, as 
suggested by the creation of vertices, £o must shrink as a grows; otherwise the discrete size 
a£ would be blown up macroscopic. Here, the evolutionary aspects capture expected prop- 
erties of full physical states, solving a vertex-creating Hamiltonian constraint, arranged by 
volume eigenvalues as internal time. (From a more general p erspec t ive, dynamical l attice s 
and their role i n quantum gravity have been discussed by IWeissI (119851 ); lUnruhl (119971 ); 
Jacobson ( 2000l ).) Then, also 5(h) must not be constant but rather a decreasing function. 



(Just as in (IT5"1) . for the regularization of the isotropic F£ b , or the c 2 in the Friedmann 
equation, by <5~ 2 sin 2 (<5c), we use the coordinate area t\ of loops, not the geometrical area 
£^a 2 which might well be constant as a changes.) 

In summary, lattice refinement means that 5 depends on the size of the universe, and 
5(h) is not a constant. T he resulting diff e rence e quations, as most generally formulated 
for anisotropic models by IBojowald et al.l ( 12007d ). are not equidistant and do not allow 
simple restrictions to separable sectors. A large class of different dynamics can thus be 
formulated in each single modelj^l Consistency conditions do exist restricting the freedom 
eve n in homogeneou s models, but systematic investigations have only just begun, as e.g. 
by IBojowald! (120091 ) . While specific details of the difference equations of loop quantum 
cosmology are not yet fully determined, owing to quantization ambiguities, there are several 
key generic features implying further properties. They are brought out most clearly in 
sovable models. 



9 Several improvised models have been formulated and sometimes analyzed in detail, based on further 
ad-hoc assumptions to constrain ambiguities. Since such an approach cannot capture the generic behavior, 
specific details of the results may not be reliable. 
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3.2.4 Harmonic cosmology 



Sometimes, solutions to the difference equations of loop quantum cosmology, which are 
linear but with non-constant coefficients, can be found. They may be studied numerically, 
and lead to int eresting issues especially in ca s es where they are not equidistan t as investi- 



gated so far by lSabharwal and Khannal (120081 ) ; iNelson and Sakellariadoul (120081 ). But in all 



cases, computing observables from the resulting wave functions and arriving at sufficiently 
generic conclusions is challenging. 

An effective treatment at this stage becomes much more powerful. Initially, one may 
expect technical difficulties due to the non-linear and non-polynomial nature of the Hamil- 
tonian constraint obtained after the loop replacement of c 2 . Fortunately, isotropic loop 
quantum cosmology offers an exactly solvable system where, in a specific factor ordering 
of the constraint operator, the dynamics is essentially free. From the point of view of 
quantum dynamics, this model is closely related to the harmonic oscillator of quantum 
mechanics, which forms the basis of most of the perturbation-theory framework for inter- 
acting quantum field theories. Similarly, the harmonic model of loop quantum cosmology 
plays a crucial role for pert urbations in quant um gravity. 

Specifically, according to iBojowald (l2007bl Jah. we obtain the solvable model for spatially 



flat isotropic loop quantum cosmology with a free, massless scalar, whose loop-modified 
Hamiltonian constraint is 



4ttG 



p 3/2 



VP1~ 2 8(P)~ 2 sin 2 (6(p)c) ~ ^c 2 + 0(c 4 ) . 



Here, S(p) (or ftf(p) = S(p)~ 3 ) is allowed to depend on the triad, incorporating different 
refinement schemes of the discrete structure. For power laws 5(p) oc p x , we introduce 
V = p 1 ~ x /(l — x) and J = p l ~ x exp(ip x c) as non-canonical basic variables forming an 
sl(2,R) algebra 

[V, J] = hJ , [V, .P] = -h.P , [J, J f ] = -2hV 

upon quantization. Thanks to the free scalar, the dynamics is controlled by the deparam- 
eterized Hamiltonian in 0, p^ cx ||?(J — jt)| =: H. 

This Hamiltonian is linear for all x, which for our linear algebra of non-canonical basic 
variables implies solvability. (Again, the absolute value is irrelevant for effective equations.) 
Following our general analysis of effective descriptions, we can directly jump to equations 
of motion for expectation values of the basic variables. They are coupled to each other, but 
not to quantum variables such as fluctuations as it normally occurs in interacting quantum 
systems. Instead of a non-linear set of infinitely many differential equations, we have a 
finitely coupled set of linear equations. For the basic variables, we have 

dji) ([./,£]) _ i/ f M_d(jt) 

dd> ih 2W d6 ' 
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As a difference to the previous discussion of effective systems, we now have to take 
into account the complex nature of our variable J. Reality conditions are required for 
physical results, which would be implemented by the physical inner product in a Hilbert- 
space representation. Here, we are not using states but directly work with expectation 
values and moments. The adjointness relation JJ^ = V 2 for our basic operators implies, 
upon taking an expectation value, a reality condition relating |(J)| 2 — (V^) 2 to moments 
of V and J. It turns out that the specific combination of moments involved is constant 
in ( ^-evolution, and that it is of the order h for a semiclassical state. (For more details, 
sec iBojowaldl (l2007al . l2008br ).) Just requiring that the state is semiclassical only once, for 
instance at large volume, ensures that the reality condition reads \ (J)\ 2 — (V) 2 = 0(h) at 
all times. With this condition, the general solution is 

(V)(0) = (£)cosh(0-/3) , ( J) (0) = -(#>(sinh(0 -/?)-*) 

with the conserved (H) and another integration constant 0. For a large class of states 
we thus have an exact realization of a bounce, with the volume bounded away f rom z ero. 
(Effective equations thus easily confirm the numerical results of lAshtekar et al.l (120061 ). at 
least as far as physical expectation values ar e conce r ned. F or the dynamical behavior of 
moments, the generic treatment of lBoiowaldl fl2007cl . l2008bl ) based on effective equations 
shows crucial differences to the numerics done for only a specific class of states.) 

However, the system is harmonic, and so its behavior is not easily generalizable to 
realistic models. While it is interesting that bounce models can be derived in this way, 
taking this at face value without realizing the solvable nature of the model would provide 
a view with severe limitations. After all, the harmonic oscillator does not provide general 
insights into quantum dynamics, and free quantum field theories allow no glimpse at the 
rich features of interacting ones. Only a systematic analysis of perturbations around the 
free model, starting with interaction terms and then introducing inhomogeneities, can 
provide a clear picture. 



3.2.5 Quantum Friedmann equation 

As already briefly discussed in the context of Wheeler-DeWitt quantum cosmology, effec- 
tive equations of quantum cosmology are necessarily of higher dimension than the classical 
equations. Quantum degrees of freedom such as fluctuations couple to expectation val- 
ues in non-solvable models, and no adiabatic or other regime has been determined where 
quantum variables could completely be expressed in terms of effective contributions de- 
pending only on the expectation values. Quantum variables are true degrees of freedom, of 
significance for the dynamics. Without a clear vacuum state to expand around, moreover, 
suitable states are less restricted than they are for low-energy effective actions. Via initial 
values for moment equations, this leaves the state dependence as a crucial contribution of 
equations, to be taken into account for sufficiently general conclusions. Since not much 
is known about details of the quantum state of the universe, results should be sufficiently 
insensitive to its properties. 
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Perturbation theory by the moments, including them order by order, is required for any 
model which is not harmonic. Examples are models with non-vanishin g spatial curvature 
or a cosmological constant a s anal yzed by iBojowald and Tavakoll ( 120081 ) (or numerically by 
Bentivegna and Pawlowskil ( 120081 )). Another class of important models is that where the 
scalar is no longer free or at least acquires a mass term. We are then dealing with a time- 
dependent Hamiltonian in the 0-evolution, and </> may not even provide a good internal 
time if the potential leads to turning points of <p{r). At this stage, a more general analysis 



i s requi red which is now available in the form of effective constraints as per lBojowald et al. 



( 12009ai ) and lBojowald and Tsobanjanl (120091 ): we do not have to deparameterize the system 
before quantizing it or computing effective equations. We can directly compute effective 
constraints and analyze them to find the physical quantum phase space. As one of the 
results, we can apply the deparameterized framework even if there is a potential provided 
it changes sufficiently slowly. This is exactly the case of interest for inflationary or other 
early-universe cosmology, and so we can derive effective equations for these situations. 

In the presence of a potential, the Friedmann equation receives the following quantum 
corrections: 
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Pcrit / ^ V P crit P 

as derived by Bojowald ( 2008dl lbl). (The first ter m, corresponding only t o the so lvable 
model with a free scalar, was earlier obtained by ISingh and Vanderslootl (120051 ): ISingh 
(|2006b. 't sums up a ll high er-order corrections due to holonomies, as explicitly expanded by 
Banerjee and Date! (120051 ). In the usual terminology of quantum field theory, it represents 
the tree- level approximation since its main form is independent of quantum back- reaction.) 
Here, P is pressure, rj parameterizes quantum correlations and 



PQ == P + 



eoPcrit + (p - P) j^j e k+1 (j^^j 



is a quantum-corrected energy density with fluctuation parameters e&. The critical density 
Pcrit = 3/87rG7 2 <5 2 V^ 3 a 2 results from the loop quantization, bringing in the scale 5(a) 
whose form is determined by lattice refinement. 

This equation contains all quantum variables in 77 and pq, subject to their own dynam- 
ics. Understanding the behavior of generic universe models in loop quantum cosmology 
requires a high-dimensional dynamical system to be analyzed. Only a few simple cases 
allow strict conclusions about the presence of a bounce. If p = P (a free, massless scalar) 
we recover the solvable scenario. For p + P large (large p^, i.e. kinetic domination), all 
corrections involving quantum variables are supressed; at least for a certain amount of 
time, the behavior can be expected to be close to the solvable one. 

In general, however, quantum variables may yank the universe away from the simple 
bouncing behavior of the solvable model. While the funda mental dynamics based on dif- 
ference equations for states remains non-singular, following IBojowald! (j2001at ). the general 
effective one is still unclear. Volume expectation values could asymptotically approach a 
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non-zero constant without bou ncing back to large values, a picture which m ight resemble 
that of the emergent universe of Ellis and Maartensl (120041 ); lEllis et al.l (120041 ) . Most likely, 
generic evolution would bring us to a highly quantum regime where effective equations or 
simple near-solvable models break down. Ultimately, the fate of a universe can be studied 
only based on the fundamental difference equations, but effective equations can show well 
how such severe states are approached. 



3.2.6 The interplay of different quantum corrections 

Loop quantum cosmology leads to different types of quantum corrections in its effective 
equations: holonomy corrections, inverse triad corrections, and quantum back-reaction. 
Depending on parameters, they may play different roles in any given regime, and some of 
them might dominate. What exactly happens can be determined only with a sufficiently 
general parameterization of correction terms. Here, all ambiguities and effects such as 
lattice refinement must be considered. 

Especially holonomy and inverse triad corrections are often closely related to each other: 
they both result from the discrete geometry, although in different ways. Comparing them 
can thus provide more restrictions on free parameters of the full theory than any one of 
the corrections would allow individually. They are both related to the size of elementary 
building blocks of a discrete geometry. For a nearly isotropic distribution in a region of 
coordinate size Vo, the classical volume on the left-hand side of 

is replaced by the right-hand side in the discrete picture. This elementary relationship first 
tells us that, dynamically, there are two free functions in the discrete picture for one free 
function a((f>) in the continuum picture: the size v(4>) and number N{<p) of discrete sites. 
Both of them typically change in internal time: the lattice they describe is being refined 
as time goes on. 

As derived earlier, one of these parameters, v, enters the basic expressions for holonomies 
exp(^ c/U 1/3 ) = exp(ic/Af 1/3 ) = exp(i 7 t; 1/3 a/a) 

l/3~ 1/3 

using c — V c = V ja. Inverse-triad corrections in isotropic models provide a correc- 
tion factor a depending on p/J\f 2 ^ 3 = (Vq/ M) 2 ^ 3 a 2 = f 2//3 . For v 1 ^ 3 ~ £p, inverse-triad 
corrections differ strongly from the classically expected value a = 1 (while holonomies 
remain nearly equal to the classically linear a /a for Hubble distances much larger than 
v 1 / 3 ). Here, quantum corrections can often be constrained. Since two different corrections 
depend on one parameter, their interplay can provide important synergistic effects in ruling 
out possible values, or even entire corrections^ 

10 We can also note that all quantum geometry corrections depend only on» = Vqo? /N ' ■ As the size 
of discrete building blocks, it is insensitive to changing Vo (both Vo and Af change proportionally to Vo) 
or coordinates (Vba 3 and TV are scaling invariant). Just like the classical equations, quantum corrections 
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Phenomenologically, v shows up in the critical density of the effective Friedmann equa- 
tion: p crit = 3/87rG7 2 f 2 / 3 . At this value, we would find the bounce of a universe sourced 
by a free, massless scalar, and even in other models its value affects details of the dynam- 
ics. This density is often assumed Plancki an, which mean s that v ~ lp. (From black- 



hole e n tropy calculations, as develo p ed by lAshtekar et al.l (119981) . iKaul and Majumdar 



(119981 ) . iDomagala and Lewandowskil (120041 ) and iMeissnerl (120041 ) . 7 ~ 0.2.) But then, 



inverse triad corrections are large. If v is constant, which is also often assumed, inverse- 
triad corrections would be large at all times, even at large total volume. A Planckian 
constant value of v is thus clearly ruled out, not by holonomy corrections alone but by 
a consistent com bination with inverse-triad effect s. Further phenomenolog i cal analysis 
i s in p r ogress, by Nelson and Sakellariadou ( 2007a bh; Copeland et al. ( 2009 ); Mielczarek 
(120081 ): iBarrau and Grain! ( l2009h : iGrain et all ( l2009bl . l2010h : IShimano and Haradal fj2009f T 
Much stronger consistency conditions can be expected when inhomogeneities are included, 
to which we will turn next. 



4 Consistent dynamics 

A space-time covariant theory is a gauge theory with the gauge group given by space- 
time diffeomorphisms. As a consequence, the local conservation law V^iG^ — 8tiGT^ v ) 
holds for the Einstein tensor G^ v and the stress-energy tensor T^ v of matter. Not all 
components of Einstein's equation are thus independent, and some can be derived from 
the others. For consistency, effective equations of quantum gravity must show the same 
form of dependencies, or allow the same number of local conservation laws, or else degrees of 
freedom would be too much constrained to have a chance of providing the correct classical 
limit. 

Consistency is automatically satisfied if one has space-time covariant higher-order cor- 
rections in Lagrangian form. At the Hamiltonian level, however, covariance is more difficult 
to check due to the lack of a direct action of space-time diffeomorphisms. There is rather 
a splitting into spatial diffeomorphisms acting as usual, and an independent generator 
deforming spatial constant-time slices within space-time. As in relativistic systems, all 
generators are constraints, taking center stage in a canonical analysis: The conservation 
law implies that the components Gq u — 8ttGT 0l/ of Einstein's equation are only of first 
order in time derivatives since <9 (67° — 8nGT®) must equal terms of at most second order 
in time. The same equation shows that these constraints are automatically preserved in 
time, given that <9o(G7° — 8ttGT®) depends only on the Einstein equations at fixed time and 
their spatial derivatives. If Einstein's equation holds at one time, time derivatives of the 
components Gq u — 8tcGT 0u vanish; the constraints can consistently be imposed at all times. 
Also this condition must be preserved for effective equations. Otherwise, the system is 

are invariant under these rescalings. For the scale factor, inverse-triad corrections become significant at a 
characteristic scale a* = (A/'/Vb) 1 / 3 ^p related to the Planck length. While a* is not scaling independent, 
it scales in the same way as a. Comparisons such as o > o, or s < a, to demarcate the classical and the 
strongly quantum regime are thus meaningful. 



29 



overdetermined with more equations than unknowns and not enough consistent solutions 
would result. 

From this perspective, the covariance condition takes the following form at the Hamil- 
tonian level. For a generally covariant system, evolution in coordinate time is generated by 
the constraints, and the constraints must themselves be preserved in time. Thus, the gauge 
transformation 5^H[N V ] = {H[N%H[e^]} of a combination H[N V ] := J d 3 xN v (G 0v - 
8ttGTq u ) must vanish for all space-time vector fields N u , e M when the constraints are sat- 
isfied. These Poisson-bracket relations provide algebraic conditions for the constraints to 
be consistent, forming a so-called first-class system. Specifically for general relativity, we 
have 

{H[N% H[M U ]} = H[K»] (21) 

with K° = C M *N° - C N aM° and K a = C M *N a - C N aM a + q ab (N°d b M° - M°d b N°). It 
is a concise expression for space-time covariance, and is insensitive to the specific dynam- 
ics: the algebra is the same for all higher-curvature actions (or with all kinds of matter), 
even though the specific constraint functionals do cha nge. (A recent Hami ltonian analy- 
sis of higher-curvature actions has been performed by iDeruelle et al.l (120091 ). showing the 
constraint algebra explicitly.) 

The canonical consistency requirement of an anomaly-free, covariant theory then states 
that the algebra of effective constraints must remain first class. As long as this algebraic 
condition is satisfied, the theory is completely consistent. The system of constraints must 
remain first class after including quantum corrections, but the specific algebra may change. 
Canonically, the realm of consistent theories is larger than at the Lagrangian level: While 
corrections to the action may be of higher-curvature form, they all give rise to exactly the 
same constraint algebra. The Hamiltonian level, on the other hand, allows changes in the 
algebra as long as it remains first class; the theory may be consistently deformed. Thereby, 
changes to the quantum space-time structure can be captured. Whether or not non-trivial 
consistent deformations exist, and how their covariance can be interpreted, is a matter of 
analysis in specific quantum gravity models. We will illustrate the results of several models 
in loop quantum cosmology in the following section. 



4.1 Cosmological perturbations 

The general issues of constrained systems in gravity can easily be seen already for linear 
cosmological perturbations. For scalar modes, perturbing the lapse function as N(t)(l + 
<f)(x,t)) and the scale factor as a(t)(l — ip(x,t)), and with a scalar matter source, the 
linearized components of Einstein's equation, in conformal time and longitudinal gauge, 
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read 



d c (ip + H(f?) = 4wGipd c 5v (22) 

VV-3-H(V + -H0) = 47rG(05^-0 2 (j) + a 2 V ip (ip)5^ (23) 

*P + H (2^ + 0) + {2H + H 2 ^ = 4nrG (<p5<j> - a 2 V tV (<p)5(p) (24) 
<9 C <9 C (0-^) = (25) 

with the conformal Hubble parameter "H. This set of equations is accompanied by the 
linearized Klein-Gordon equation for tp — (p + 5<p. The first two lines are of first order in 
time, and thus pose constraints for initial values. They are the linearized diffeomorphism 
constraint ( )22|) and the Hamiltonian constraint ( 1231) . The next two lines are the diagonal 
and off-diagonal components of the spatial part of Einstein's equation, also linearized 
around Friedmann-Robertson-Walker. Using the background equations for % and the 
background scalar (p, one can explicitly derive the Klein-Gordon equation from the rest. 
The system is thus overdetermined, but in a consistent way. Canonically, this follows from 
the fact that all equations result from a set of first-class constraints. 

In these equations, perturbations 0, ip and 5<p are directly used for the metric and 
the scalar. These quantities are subject to gauge transformations under changes of space- 
time coordinates, which also change the gauge. So far, the equations as written are in 
longitudinal gauge where the space-time metric remains diagonal under perturbations. 
More generally, scalar modes can be perturbed also in the off-diagonal part of the spatial 
metric by Sq a b = d a dbE, or in the space-time part by a perturbed shift vector 5N a = d a B. 
Two new scalar perturbations E and B are introduced in addition to and if), and all 
four transform into each other under coordinate changes. For linear coordinate changes, 
however, the combinations 

* = -0 - n(B - E) (26) 

and 

$ = 0+ (b-e)' + H(B-E) (27) 



and a similar one for the matter perturbation remain unchanged as specified by iBardeen 



(119801 ); they form gauge-invariant observables of the linear theory. From a canonical per- 
spective, these combinations are invariant under the flow generated by the constraints. 

The ungauged evolution equations, containing E and B in addition to and ip, can 
be expressed completely in terms of the gauge-invariant variables, as required on physical 
grounds. Here, the first-class nature of constraints is important which ensures that the 
constraint equations, and thus ff22|) and ff23l . are gauge invariant. Consistency as well 
as gauge invariance of the equations of motion is thus guaranteed by having a first-class 
algebra of constraints. Also here, the first-class nature serves as a complete requirement 
for covariance. 
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4.2 Gauges and frames 

A first-class algebra of constraints ensures that physical evolution can be formulated fully 
in gauge- invariant terms. By the same property of the constrained system, also frame- 
independence is guaranteed: The gauge algebra corresponds to deformations of the spatial 
hypersurfaces given by a time function t = const. Gauge invariant variables are insensitive 
to these deformations, and thus to the choice of time. But a time function (without speci- 
fying which spatial coordinates are to be held fixed) does not uniquely tell us how to take 
time derivatives for the equations of motion; this is only accomplished when we also specify 
a time-evolution vector field t a (such that t a d a t = 1). For a given time function, there are 
many choices for t a , and the time-evolution vector field may be changed independently of 
the foliation. 

A fixed foliation of space-time into spatial slizes with unit normals n a allows us to 
express the freedom in choosing time-evolution vector fields by the lapse function N and 
the shift vector N a (with N a n a = 0) such that t a = Nn a + N a . These are exactly the 
functions which appear in the constraints generating evolution: for a fixed choice of t a , 
or N and iV a , Hamiltonian equations of motion are / = C^f = {/, H[N, N a }} for any 
phase-space function /. Since N and N a appear as multipliers of first-class constraints, 
they can be chosen arbitrarily (except that we would like N > for evolution toward 
the future). For a consistent set of first-class constraints, we can thus choose the frame 
freely, and different frame choices are guaranteed to produce consistent results. (In a 
space-time treat ment, observables which are gauge as well a s frame-independent can be 



derived following Ellis and Brunil (119891 ): iBruni et al.l (119921 ) . In a reduced phase-space 



treatment, where one uses observables solving the classical constraints, one is working at 
a gauge-invariant, but not automatically frame- independent level.) 

In classical relativity, cosmological perturbation equations can often be derived in much 
simper ways when a space-time gauge is chosen, such as the longitudinal gauge above or 
the uniform one where only matter fields are perturbed. Since gauge transformations are 
known to correspond to space-time coordinate transformations, one can directly verify that 
such gauges are possible. Moreover, choosing a gauge before deriving equations of motion 
from an action or Hamiltonian is equivalent to choosing a gauge in the general equations 
of motion. The situation is thus completely unambiguous. 

When quantum effects are included, either in a full quantum theory or in an effective 
manner, the constraints change by quantum correction terms. Equations of motion change, 
as expected, and so do the form of gauge-invariant expressions since it is the constraints 
which generate gauge transformations. In such a situation, gauge transformations and 
suitable gauge fixings can be analyzed only after the quantization has been performed and 
the corrected constraints are known. If gauge fixings are employed before quantization or 
before determining effective constraints, the choice of gauge fixing may not be compatible 
with the resulting corrected gauge transformations. Moreover, choosing different gauge 
fixings before doing the same kind of quantization would in general lead to different final 
results, making the procedure ambiguous even beyond unavoidable quantization ambigui- 
ties. Similarly, a reduced phase space quantization is based on frame-fixing, although no 
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gauge need be fixed. 

For the different approaches u sed in relation to loop quantum gravity, several exam- 



ples exist. ICampiglia et al.l (120061 ) develop methods to deal with a discretization possibly 



breaking gauge symmetries. Similar methods have then be used by lCampiglia et al.l (120071 ) 
in spherically symm etric models with (pa r tial) g auge fixing. Also the hybrid quantizations 
of Gowdy models bvlMartm-Benito et al.l (120081 ) rely on gauge fixing of the inhomogeneous 



generators. iBahr and Dittrichl (l2009bl ) const ruct discretized theories in three dimensions, 



respecting the space-time gauge, but argue i nlBahr and Dittrichl (l2009al) that this may not 



be po ssible in four space-time dimensions. lLaddhal (120071 ) and lLaddha and Varadarajan 



(120081 ) quantize 2-dimensional parameterized models of field theories by the Dirac proce- 
dure and represent observables on the resulting physical Hilbert space. In this treatment, 
the discretization does not lead to inconsistencies but possibly to deformations of clas- 
sical algebras of observables. Finally, reduced pha se space methods fix ing the frame by 
referring to an extra dust field are deyelope d by iGiesel et al. (l2007al Jbr) for cosmologi- 
cal perturbations, and by IGiesel et al.l (120091 ) i n spherical symme t ry. O ther treatments 
for_cosmo logical perturbations are u sed, e.g. , by Artymowski et al.l ( 2009 ) and Mielczarekl 
(120101) in gauge-fixed vers i ons an d by lPuchtal (120091 ) in a frame- fixed (reduced phase-space) 
way. Rovelli and Vidotto ( 20081 ) provide a proposal by which cosmological perturbations, 
taking into account space-time discreteness, might be implementable by a consistent first- 
class algebra of constraints. So far, this has been realized for coupling two independent 
homogeneous patches. 

The only valid treatment of a complicated gauge theory is by working without restric- 
tions of the gauge throughout the quantization procedure, until the final gauge algebra 
resulting from the corrected constraints has been confirmed to be consistent. Here, the 
anomaly problem, confirming that a consistent deformation is realized, must be faced head- 
on and cannot be evaded. In the final equations one may choose one of the allowed gauges 
for further analysis, but gauges cannot be used to simplify the quantization. In the rest of 
this exposition, we follow these lines to illustrate the consistency of several effective sets of 
constraints incorporating some of the discreteness effects of loop quantum gravity. 



5 Consistent effective discrete dynamics 

If loop quantum gravity has a chance of being a viable quantum theory of gravity, the form 
of discrete quantum geometry it implies must give rise to effective dynamical equations 
satisfying the consistency conditions of a covariant theory. At the Hamiltonian level, this 
requires a first-class algebra of constraints. After the preparation in the preceding sections, 
we can now see what specific models indicate. 

5.1 Constraint algebra 

Consistent deformations implementing the effects of loop quantum gravity have been found 
in several different cases. Most of them use inverse-triad corrections, which have been in- 
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corporated successfully in spherically symmetric models by iBojowald and Reyesl (120091 ); 
Bojowald et al.l (l2009d ) as we ll as linear perturbati o ns aro und spatially flat Friedmann- 
Robert son- Walker models by IBojowald et al.l (12008L l2009bl ). The situation for holonomy 
corrections is more restric t ive; h ere, certain versions have been realized in spherically sym- 
metric models by iReyed (120091 ) as well as for linear tensor and vec t or mo des around 
Friedmann-Robertson- Walker models in IBojowald and Hossainl (120071 . 120081 ) . However, 
so far no inhomogeneous model has been found where holonomy corrections in a com- 
plete form would be consistent. Here, the requirement of anomaly- freedom appears very 
restrictive. 

Inverse-triad corrections have been implemented consistently in several settings, and 
so are not ruled out by consistency considerations. In particular, corrections from the 
discreteness of quantum geometry are allowed and do not necessarily spoil covariance. The 
specific form of their implementation then tells us if and how space-time structures have 
to change due to quantum effects. 

Based on formulas such as (Il4p . inverse-triad corrections arise for any term in the Hamil- 
tonian constraint bear ing components of the inverse densitized triad, such as 1/ \J | det Ef \ . 
Bojowald et al.l (120081 ) have consistently implemented these corrections for linear inhomo- 
geneities, where the corrected constraint algebra is of the form (121]) but with 

K a = £ M aN a - C N aM a + a 2 Na- 1/2 d a {5M° - 8N°) 

while K° = £m*8N — C^aSM retains its classical form. In addition to the contribution 
Na~ 1 / 2 d a (SM° — 8N°), which is expected classically for a linearization around Friedmann- 
Robertson- Walker models with iVj = N + 8Ni, there is the function a (depending on the 
background scale factor a) arising from inverse-triad corrections. An algebra of the same 
for m aris e s for spherically symmetric models, with different matter couplings as discussed 
by iReyesI (120091 ). The constraint algebra is anomaly-free: the system of constraints remains 
first class. But it is not exactly the classical algebra, and thus deformed. Inverse-triad cor- 
rections from loop quantum gravity cannot amount to higher-curvature corrections to the 
action since this would leave the constraint algebra unchanged. Rather, these corrections 
can only be understood as deforming local space-time symmetries while keeping covariance 
realized. 



5.2 Cosmological perturbations 

Consistent versions for quantum-corrected constraints allow one to analyze their impli- 
cations for the dynamics. When constraints are corrected, not just evolution equations 
change but also the gauge transformations generated by the constraints. Thus, expressions 
for gauge-invariant observables depend differently on perturbations of the fields, which by 
itself may give rise to new effects. Other implications then follow from studying the dy- 
namical evolution of gauge-invariant observables. (Some quantum-gravity corrections have 
been implemented in gauge-fixed versions. They are formally consistent, but quite arbi- 
trary in their implementation. For instance, it remains unclear how different gauge-fixings, 
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all done before quantization, might affect the results. Moreover, some physical effects due 
to corrections to gauge-invariant observables can easily be overlooked.) 

For linear perturbation equations around Friedmann-Robertson-Walker models, cos- 
mological perturbation equations are the main application of consistent deformations. With 
a consistent deformation, perturbation equations form a closed set and can be written 
fully in terms of gauge -invariant variables. For inverse-triad corrections, as derived by 



Bojowald et al.l (l2009bf ). they take the form 

d c + n{\ + /)$) = icG^4>d c 6cp 

as the corrected time-space part of Einstein's equation, 

A(a 2 ^) - mil +f)(if + HQ(1 + f) 



= 4ttG|(1 + / 3 ) (^ GI - f(l + A)$ + va 2 V v {(p)5 V Gl ) 
as the corrected time-time part, and 

v + nfa (i-±f) + * { i + f)) + ( 2 n + n* + 



2a da J J \ \ 2 da 2a da J 

= 4ttG| (05if Gl - a 2 vV v (£)5ip G1 ) 

as the diagonal space-space part. All corrections /, f±, and h below are determined from 
the basic inverse-triad corrections a (for the gravitational part of the constraint) and v (for 
the kinetic term of the matter part). Since these are background corrections, their form 
can easily be computed in isotropic models, suitably parameterized for all ambiguities and 
lattice refinement. 

The off-diagonal space-space part also implies a non-trivial equation, with an unex- 
pected consequence: While the classical analog would simply identify $ = the corrected 
equation implies $ = + h) with a quantum correction by h ^ 0. This may be inter- 
preted as an effective anisotropic stress contribution, but it results from a correction to 
quantum gravity, not to matter. 

As a secon d implication, w e may have no n-conservation of power on large scales as 



pointed out by iBojowald et al.l ( I2007bl . l2009bl ) J I This may be important for inflationary 



structure formation, where the long evolution times while modes are outside the Hubble 
radius would make even a weakly changing size of the overall power significant. Both 
of these effects are difficult to see in gauge-fixed treatments, such as the longitudinal or 
uniform gauge, or in frame-fixed versions based on reduced phase space quantizations. 



ll Note added: Recently it turned out that power, due to an unexpected cancellation, is con- 
served on large scales for the currently available equations with inverse-triad corrections, as derived by 
Boiowald and Calcagn i (2010). So far, however, there is no general argument for th e conse rvation of power 



in the presence of modified space-time structures. The analysis of iBoiowald et al. I d2010h has shown that 



interesting effects for potential observations mainly arise from corrections to the running of indices. 
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Since scalar cosmological perturbations have been consistently formulated only for inverse- 
triad corrections, no version is yet able to evolve perturbations through a bounce for 
which holonomy corrections are required. (Implementing holonomy corrections only for 
the background leads to inconsistent evolution equations for inhomogeneities.) 

In addition to inverse-triad corrections, holonomies and quantum back-reaction must 
be implemented to obtain a full picture from loop quantum gravity. While consistent 



deform ations are not yet known for the latter two corrections, iBojowald and Skirzewski 



( 120081 ) have formulated quantum back-reaction in a cosmological setting. In general, in 
quantum gravity this requires the inclusion of moments between all degrees of freedom of 
gravity and matter, including quantum correlations between them. By setting the quantum 
variables of gravity as well as its quantum correlations to zero, one obtains the effective 
equations of quantum field theory on a curved space-time as a limit. Including leading order 
corrections from the gravitational quantum variables provides quantum field theory on a 
quantum space-time. While such limiting cases can be realized explicitly at the effective 
level, much still remains to be done for a detailed analysis of specific properties. 

5.3 Causality 

Several examples illustrate the importance of a consistent constraint algebra, rather than 
just any deformation as allowed in gauge- fixed or frame-fixed treatments. We have already 
seen that some effects in cosmological perturbation equations can be obtained only when 
neither gauge nor frame are fixed before the theory is quanti zed or effective constraints are 



derived. Another example is the realization of causality, as IBojowald and Hossainl (120081 ) 
studied it by comparing the propagation of gravitational waves to that of light. 

With quantum-gravity corrections, the gravitational as well as the Maxwell Hamiltonian 
change compared to the classical expressions, affecting evolution equations and their plane- 
wave solutions. The gravitational contribution to the Hamiltonian constraint for inverse- 
triad corrections is 

H - = Tito L A3xa{E;) 70T\ ( e>itFL - 2(1 + -< 2)K '^ 

implying the linearized wave equation 

1 / 1^ + 2 a / _ . _ \ = 

2 \a a \ a J J 

for the tensor mode h l a on a cosmological background with scale factor a and source-term 
IP a . By a plane-wave ansatz, we derive the dispersion relation u 2 = a 2 k 2 for gravitational 
waves. Since a > 1 for perturbative corrections, there is a danger of gravitational waves 
being super-luminal. 

With these corrections, gravitational waves are faster than light on a classical back- 
ground. For a meaningful comparison, however, we should use the speed of gravitational 
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waves in relation to the physical speed of light on the same background, which should 
receive quantum corrections, too. Here, the Hamiltonian is 

Hem = [ d 3 x(a EM (q cd )^=E a E b q ab + (3 EM (q cd )^^F ab F^q ac q bd ) 
V Vdetg 16tt J 

with two correction functions «em and /?em kinematically independent of the gravitational 
correction a. From the wave equation 

d t (a^d t A a ) - f3 EM V 2 A a = 

we obtain the dispersion relation u 2 = ckem/Sem^ 2 , which also is "super-luminal" compared 
to the classical speed of light. 

Working out the requirements for anomaly-freedom with these two contributions to the 
Hamiltonian constraint, as done by iBojowald and Hossainl (120081 ). we find a 2 = oem/^em 
and the dispersion relations are equal. Physically, for comparisons of speeds on the same 
background, there is no super-luminal propagation. In a gauge-fixed or frame-fixed treat- 
ment, by comparison, one could have chosen the correction functions independently of 
another since tensor modes of the gravitational field and the electric field make up in- 
dependent physical observables. A gauge- or frame-fixed treatment could easily produce 
corrected equations violating causality, but this is ruled out by a complete treatment. 



6 Outlook: Future dynamics 

To probe a quantum theory of gravity or even arrive at predictions one must evaluate its 
dynamics in detail. For low-energy effects, leading corrections to classical equations must 
be derived. The best tool for systematic investigations in such cases is that of effective 
descriptions, providing the evolution of expectation values of observables in a physical 
state. At the same time, they can tell us much of the entire behavior of physical states. 

Hamiltonian effective descriptions can be applied directly to canonical quantum grav- 
ity and cosmological models. In particular, typical implications such as the discreteness of 
spatial or space-time structures can then be probed, or first ensured to be consistent at all. 
Several examples in loop quantum cosmology have demonstrated that discreteness correc- 
tions can indeed be implemented consistently, leaving the theory covariant but deforming 
its local space-time symmetries. Future work must ensure that this is indeed possible for 
the full theory of loop quantum gravity and its effective constraints. 

While isotropic solvable models of loop quantum cosmology suggest a role of bouncing 
cosmologies for potential scenarios, no consistent set of equations to evolve inhomogeneities 
through a bounce has been found. The only available options so far make use of gauge 
(or frame) fixings before quantization, and thus miss crucial aspects of space-time struc- 
tures. Any mismatch of growing modes in the collapse and expansion phases can easily be 
enhanced by cosmic evolution, providing opportunities for potential observations but also 
requiring extreme care in finding fully consistent equations. Inhomogeneous cosmological 
scenarios remain uncertain, and with it follow-up issues such as the entropy problem. 
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What models investigated so far suggest, in many different versions, is that the clas- 
sical algebra of space-time diffeomorphisms is deformed but not violated. It is then clear 
that quantum corrections cannot merely amount to higher-curvature terms in an effective 
action, although such terms may appear, too. Instead, quantum structures of space-time 
must change by quantum effects. While space-time covariance is no longer realized in the 
standard sense, from the Hamiltonian perspective the effective theories remain completely 
consistent and covariant with an underly i ng fir st-class algebra of gauge generators. As 
indicated by the initial quote from iDirad (119581 ) . this sense of covariance, from a general 
perspective, is the appropriate one. Its consistent implementation, without fixing gauge or 
frame, can tell us a great deal about the fundamental structures of space and time. 
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